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ABSTRACT 


Chemistry educators teaching chapters on "real gases" usually present to the students the "virial development" of Kamerlingh Onnes for the equation of state, which 
takes into account interactions of the molecules in aggregates by a successive series of parameters (B, C, etc.), yielding a polynomial equation on the concentration 
(equivalent to the inverse of the molar volume). A well known transformation of the virial equation to a polynomial equation as a function of pressure is also widely 
used, with suitable coefficients (B', C', etc.). The literature presents relationships between the second and third virial coefficients for both equations, for example B' = 
B/RT. Students usually assume that both formulations should deliver to similar calculated results. However, when using the relationships mentioned above to perform 
calculations for the same gas conditions of pressure and temperature, discrepancies are observed in the molar volume obtained by the two formulations of the virial 
equation. Chemistry educators are required to present such discrepancies, explain them and interpret them. In addition to that, a possible alternative practical approach 
is to develop a relationship between the second virial coefficients B and B', that yields suitable results when using both formulations. 
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The "virial expansion" is considered one of the most interesting and versatile 
equations of state for gases(Salzman, 2004). It was developed by Kamerlingh 
Onnes at the end of the 1 9th century as an improvement of the van der Waals equa- 
tion, based on the law of corresponding states(Wisniak, 2003), and it is used 
almost in every Physical Chemistry course as one of the possible tools describing 
behavior of "real" gases.The virial expansion is a power series in powers of con- 
centration, and has the general form: 


P = cRT(l+ 3c + Cc 2 + 



( 1 ) 


where P is the pressure, c the molar concentration, R the gas constant, T the abso- 
lute temperature, B and C the second and third virial coefficients, respectively, 
and V m the molar volume, defined as the reciprocal of the concentration, thus V m 
= 1/c = V/n, where V is the volume and n the amount of moles. The successive 
virial coefficients B, C, etc. take into account deviations from the ideal gas law 
due to interaction of the molecules in pairs, triplets, etc(Riddell & Uhlenbeck, 
1953). 

The virial equation truncated at the second coefficient is considered a relatively 
simple, yet accurate method for describing the state of gases, even at relatively 
high pressures(Hayden & O’Connel, 1 975). Considering several presented theo- 
retical and empirical methods to evaluate the second virial coefficient(Dymond, 
Cholinski, Sz^franski, & Wyrzykowska-Stankiewicz, 1986; Hayden & 
O’Connel, 1975; O’Connel & Prausnitz, 1967; Pitzer & Curl, 1957; Riddell & 
Uhlenbeck, 1953; Tsonopoulos, 1 974), a useful equation of state for low or mod- 
erate vapor densities may be based on determination of the required B parameter: 
when the temperature and the molar volume are known, the pressure might be 
evaluated using: 


B' and C' are also known as virial coefficients. The relationship between these 
parameters and the "original" virial coefficients is usually denoted as(Atkins & 
de Paula, 2006; Boschi-Filho & Buthers, 1997; Salzman, 2004): 



, C r = 


C-B 2 

{rt) 2 


( 5 ) 


Students usually will assume that by using the truncated virial equation (Eqs. [2], 
[4]), as is commonly done (Hayden & OConnel, 1975; OConnel & Prausnitz, 
1967; Pitzer & Curl, 1957), combined with the relationships presented in Eq. [5] 
should deliver the same P, V and T values. However, it is the duty of the Chemis- 
try Educators to point out that in fact this practice yields inconsistencies in the 
evaluated values. For example, for C0 2 at 280K, B = -0.146 x 10' 3 m 3 mol' 1 
(Kehiahian, 2005). Evaluating the molar volume of C0 2 at 280K with a pressure 
of 2 MPa by transforming Eq. (2) into a quadratic equation and solving for V m 
yields 0. 9926 x 10' 3 m 3 mol' 1 , andZ= 0.853. On the other hand, calculation of B' 
= -6.354 x 10' 8 Pa' 1 using Eq. [5] followed by Eq. (4) yields V m = 1.018 x 10' 3 m 3 
mol' 1 and Z= 0.875 — a discrepancy of about 3%. Discrepancies increase consid- 
erably with pressure, and at 3.5 MPa, differences exceed 15%, with Z values of 
0.675 and 0.780 for Eqs. (2) and (4), respectively. 


The use of the less truncated equation, including the third virial coefficient con- 
siderably improves the fit between the two formulations, but even then, discrep- 
ancies remain: for example, very detailed values of B and C are given by Estrada- 
Alexanders and co-workers from the integration of speed of sound (Estrada- 
Alexanders, Guzman, & Perez- Vidal, 2012). For C0 2 at 280K they report B = - 
0. 142 1 x 10' 3 m 3 mol' 1 and C= 5.278 x 1 O' 9 m 6 mol' 2 . By using Eq (5), values of B'= 
-6.105 x 10' 8 Pa' 1 and C' = -2.754 x 10' 15 Pa' 2 are obtained. At a pressure of 3.5 
MPa, F m and Z values are (0.4586 x 10' 3 m 3 mol' 1 , 0.715) and (0.4759 x 10' 3 m 3 
mol 1 , 0.742) using Eqs. (1) or (3), respectively, a close to 4% discrepancy. 




( 2 ) 


A detailed examination focusing on the truncated equation shows that such dis- 
crepancies are indeed inherent: combination of Eqs. [2] and [4] delivers two dif- 
ferent definitions for the compressibility that should converge: 


where Z is the "compressibility" or "compression factor", which is a useful mea- 
sure of the deviation of a real gas from an ideal gas. For an ideal gas, the compres- 
sion factor is equal to 1 (Atkins & de Paula, 2006; Salzman, 2004). If the pressure 
is known, V m can also be evaluated by transforming Eq. (2) and solving a qua- 
dratic equation. However, where the molar volume is to be determined, an alter- 
nate form of the virial equation is usually presented as(Atkins & de Paula, 2006): 


p y , v 

2- = [l+ B' P + C P 2 ...) 

RT ' ’ 


( 3 ) 


which, when considering only the second virial coefficient yields: 


PV 

Z = ^ = {l + P'P) 

RT 


( 4 ) 



However, introducing the commonly used relationship between B and B' 
(Eq.(5)) yields: 

(■♦ 

( B 


PV m =RT 

presumably indicating that the gas behaves according to the perfect gas equation 
of state, with Z = 1 .00. 


RT 


V. 




( 7 ) 
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The discrepancy between the "reality" (P, V, and T of the actual gas) and the math- 
ematical models developed to describe such reality should be elaborated in front 
of the students. G.E.P. Box (Box, 1976) explains that "In applying mathematics 
to subjects such as physics or statistics we make tentative assumptions about the 
real world which we know are false but which we believe may be useful nonethe- 
less". 

For students seeking a coherent model that can yield the same calculated results 
for both Eqs. [2] and [4], an alternative to Eq. [5] can be developed, yielding a dif- 
ferent relationship between both "second virial coefficients". Reorganizing Eq. 
(6) leads to: 


-[pV„\-RT(l + B'P)- RT f 1 + y^ 

V / V mJ 

where the definition of B' as a function of V m can be readily used if the molar vol- 
ume is known. When pressure is the known value, as in the examples given 
above, Eq. (8) can be elaborated to isolate B' to £'+.zr 2 P = which can then 
be easily solved as a quadratic equation. RT 

Using this approach on the example presented above, and calculating the molar 
volume of C0 2 atP = 2 MPa and T= 280K, leads to B'= -7.353 x 10' 8 Pa' 1 , and the 
use of Eq. [4]yields V m = 0.9926 x 10' 3 m 3 mol' 1 , and Z= 0.853, exactly the same 
results as from Eq. [2]. The same happens for P = 3.5 MPa: B'= -9.298 x 10' 8 Pa' 1 
and V m = 0.4486 x 10' 3 m 3 mol' 1 , Z— 0.675 for both methods of calculation. 

Thus, by evaluating B 'using the relationship presented herein, consistent results 
of the two formulations of the virial equation are achieved. However, students 
(and educators as well) might recall that this is only a model, and "Since all mod- 
els are wrong the scientist cannot obtaina "correct" one by excessive elabora- 
tion. "(Box, 1976). 
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